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DYNAMIC  ANALYSIS  OF  THE  PIXIE  BUOY 

The  PIXIE  Buoy  has  been  chosen  as  the  initial  platform  to  be  used  for 
acquiring  airborne  magnetic,  electric,  and  electromagnetic  measurements  during 
the  preliminary  sea  trials  of  the  LINEAR  CHAIR  system  in  FY  80.  A  mooring  system 
must  be  developed  that  will  minimize  buoy  motion  response  in  both  shallow  (100  to 
200  feet)  and  deep  (2000  to  3000  feet)  water  configurations.  The  following  analy¬ 
tical  data  are  illustrative  of  a  technique  to  predict  and  evaluate  the  expected 
responses  of  the  PIXIE  Buoy  under  various  sea  conditions  and  mooring  configurations. 
These  equations  and  graphs  are  meant  as  an  aid  to  the  investigator  in  designing 
the  optimum  mooring  configuration  for  the  buoy  with  regards  to  allowable  motion 
and  tilt  as  determined  by  the  sensitivity  of  the  buoy-mounted  instruments.  Actual 
numerical  values  presented  are  only  illustrative  as  they  do  not  reflect  the  actual 
configuration  of  the  PIXIE  Buoy  as  it  will  be  used  in  the  LINEAR  CHAIR  tests.  Once 
the  actual  numerical  values  for  the  various  parameters  have  been  established,  the 
following  technique  can  be  used  to  generate  more  accurate  values  of  the  buoy 
response. 
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Section  1 


STEADY  STATE  PITCH 

1 .  BUOY  PITCH  VERSUS  SUBMERGENCE  AND  CURRENT  SPEED 

The  PIXIE  Buoy  is  depicted  in  Figure  1-1  as  a  free  body  with  the  steady 
state  forces  acting  on  it.  The  current  is  acting  on  the  buoy  from  the  right  hand 
side  and  the  forces  T^  and  T  represent  the  horizontal  and  vertical  components 
of  tension  which  the  buoy  imparts  to  the  mooring  line. 

Two  waterlines  are  shown  in  Figure  1-1.  The  lower  of  the  two,  the  free 
floating  waterline,  represents  the  case  where  the  buoy  is  freely  floating  in  still 
water  (i.e. ,  no  current)  and  no  mooring  line  attached.  It  is  located  a  distance  L 
from  the  bottom  of  the  buoy.  The  higher  of  the  two,  the  loaded  waterline, 
represents  the  case  where  the  buoy  is  moored  in  a  current.  The  parameter  H 
represents  the  distance  between  the  two  water  lines,  or  in  other  words,  the  buoy 
submergence. 

The  forces  D^  and  D^  represent  the  normal  and  tangential  components 

of  hydrodynamic  drag  due  to  the  steady  state  current.  For  a  uniform  current 

acting  in  the  region  of  the  buoy  (assumed  for  this  analysis)  the  drag  forces  will 

act  at  the  mid -point  of  the  submerged  portion  which  is  located  at  a  distance  of 

(L  +  H)/2  from  the  bottom  of  the  buoy, 
o 

The  forces  W  and  B  represent  the  buoy  weight  and  buoyancy  acting  at 

distances  of  Lw  and  from  the  bottom  of  the  buoy.  The  buoyancy  force,  B, 

represents  the  free  floating  or  unloaded  displacement.  The  increased  buoyancy, 

A  B ,  due  to  the  buoy  submergence  from  the  free  floating  condition  acts  at  a 

distance  of  L  +  H/2  from  the  bottom  of  the  buoy, 
o 

The  angle,  &,  is  the  resultant  buoy  steady  state  pitch  angle. 

A  summary  of  the  buoy  physical  parameters  are  presented  in  Table  1-1 


Figure  1-1.  PIXIE  Buoy  Freebody 
Diagram 


Table  1-1.  PIXIE  Buoy  Physical  Parameters 

Buoy  Weight  W  57,169  lbs 

Free  Floating  Buoyancy  B  57,169  lbs 

Buoy  Diameter  d  4.77  ft 

Free  Floating  Waterline  L  50  ft 

o 

C.G.  Location  L  15.53  ft 

w 

C.  B.  Location  L  24.99  ft 

B 

Summing  the  forces  in  the  x  and  y  directions  and  the  moments  about  the 
towpoint  we  get: 


E  F  =  0 

^  x 


T  =  D  cos  e  +  D_  sin  e 
x  N  T 


y  f  =  o 


T  =  B+aB-W  +  D_  cos  6  -  D„t  sine 
y  T  N 


Em  =0 


H  L  +  H 

BL  sin  9  +  AB  (L  +  —  )  sin  8  -  D  (  — -  )  -  W  L  sin  6  =  0 

d  O  Z  i »  z  w 


L  +  H 

(BL3  +  AB  (Lq  +  — )  -  w  Lw)  sin  8  =  Dn  (  2  ) 


6  =  Sin 


-1 


L  +  H 

dn  (  -“i - ■» 


H, 


blb  +  +  ?>  -  VVLw 


The  normal  hydrodynamic  force,  D^.,  is  given  by 


D  =  1/2  e  CM  Af  0/  cos  6  ) 


w 


where  V  is  the  current  speed,  V  cos  0  is  the  component  of  velocity  normal  to  the 
cylinder,  P  is  the  fluid  density,  C  is  the  normal  drag  coefficient,  and  A  is  the 

projected  frontal  area  given  by 

A  =  (L  +  H)d 
1  o 

Therefore ,  the  buoy  pitch  angle  is  given  by 

1/2  P  C  .  (L  +  H)  d  V2  cos  e2  (  ) 

_ ° _ 

BLb  +  AB(Lo+ 


the  increased  buoyancy,  aB,  is  given  by 


AB  = 


H 


For  small  angles  of  e  ,  then 

cos  e~l 
2 

cos  0=1 

and  the  expression  for  e  reduces  to 


9  =  Sin 


P  C.T  (L  +  H)2  d  V2 

NB  ° _ 

4BL_  +  P  g  -  d2K  (L  +  —■ )  -  4  WL 
B  o  2  w 


For  values  of  buoy  submergence,  H,  between  0  and  8  feet  and  for  current  speeds, 
V,  between  0.1  and  1.0  knots  the  resultant  pitch  angles,  0,  are  shown  in 
Figure  1-2.  The  values  of  the  constants  used  in  the  above  calculations  were 
assumed  to  be  as  follows: 


p 


1.2 

1.99  slug  s/ft2 
32.2  ft/sec2 


g 
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BUOY  PITCH  VERSUS  MOORING  LINE  SCOPE 


Until  now  the  steady  state  buoy  pitch  angle  has  been  described  strictly 
in  terms  of  the  buoy  as  a  free  body  by  itself  where  the  buoy  submergence,  H, 
has  been  a  variable  parameter.  It  has  been  shown  that  for  a  given  current  speed 
and  a  given  buoy  submergence  there  is  one  and  only  one  resultant  buoy  pitch  angle. 
The  buoy  submergence  parameter  is  not  an  independent  variable,  however.  The 
key  to  this  problem  lies  in  the  tension  vector  resulting  from  the  mooring  line 
constraint.  The  parameter  H  must  now  be  related  to  some  other  parameter  of  the 
overall  mooring  system. 

It  can  be  shown  that  for  a  given  current  profile,  water  depth, mooring 
line  scope  and  buoy,  the  mooring  system  will  assume  some  steady  state 
equilibrium  configuration.  That  is,  for  any  one  set  of  conditions  there  will  be 
one  and  only  one  resultant  mooring  line  spatial  geometry,  and  there  will  be  one 
and  only  one  resultant  buoy  attitude  in  terms  of  forces,  pitch  angle  and  submergence. 
Therefore,  for  a  given  current  profile  the  resultant  buoy  submergence  (and  therefore 
buoy  pitch)  is  a  function  only  of  mooring  line  scope. 

In  reality,  scope  and  current  speed  would  be  the  independent  parameters 
and  buoy  submergence  would  be  the  dependent  parameter.  However,  since  there 
is  a  one  to  one  correspondence  between  scope  and  submergence,  it  is  convenient 
for  the  purpose  of  this  analysis  to  reverse  the  relationship  and  make  submergence 
and  current  speed  the  independent  parameters  and  scope  the  dependent  parameter. 

The  approach  is  to  evaluate  the  relationships  for  T^  and  Ty  previously 
formulated  as  a  function  of  buoy  submergence  and  use  the  resultant  tension  vector 
as  an  initial  condition  to  the  computer  program  TOWLIN  to  calculate  the  mooring 
line  scope  required  to  traverse  the  depth  of  water.  The  mooring  configuration 
scenario  is  depicted  in  Figure  1-3. 

It  is  assumed  that  the  water  depth  is  2000  feet.  It  is  further  assumed 
that  the  surface  current  acts  only  in  the  region  of  the  buoy  (approximately  0-50 
feet  below  the  surface)  and  that  the  subsurface  current  speed  is  3  0  percent  of  the 
surface  current  speed . 
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The  mooring  line  is  assumed  to  have  the  following  properties: 


Material 
Construction 
Jacket 
Diameter 
Weight  in  Air 
Weight  in  Water 
Breaking  Strength 


Kevlar  29 

6x19  stranded  with  independent  core 

Polyurethane 

0.61  in. 

0.114  lbs/ft 
0.0  lbs/ft 
30,000  lbs 


The  normal  and  tangential  hydrodynamic  cable  drag  coefficients  were 
assumed  to  be 


CN  '  I*S 
c 


CT  =  .005 

c 


As  previously  formulated,  the  resultant  components  of  tension,  T  and 


Ty,  at  the  buoy  are  given  by 


T  =  D  cos  0  +  D  sin  6 
x  N  T 


T  =  p  g 
y 


— —  H  +  Dt  cos  0  -  sin  B 


where 


D„  = 


'N 


1/2  o  C„  (L  +  H)  d  V2  cos2  g 
NB  0 


DT  " 


1/2  P  C  (L  +  H)  -  d  V2  sin2  e 
TB  0 


e  = 


,  -l 
sin 


P  C 


NB  (Lq  +  H)2  d  v2 


4  BL  +  p  g  -  d2H  (L  +  —  )  -  4WL 
B  o  2 


w 
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Using  T  and  T  ,  the  initial  conditions  to  the  mooring  problems,  the 
x  y 

computer  program  TOWLIN  numerically  integrates  the  differential  equations  which 
describe  the  tension  in  the  cable  and  the  curvature  of  the  cable  due  to  the 
hydrodynamic  and  weight  forces  acting  on  an  incremental  length  of  cable.  The 
program  proceeds  with  the  integration  until  the  specified  depth  of  2000  feet  is 
reached.  The  resulting  cable  scope  as  a  function  of  buoy  submergence,  H,  for 
two  current  profiles  is  presented  in  Figure  1-4, 

It  can  be  seen  that  there  is  a  strong  dependence  of  buoy  submergence 
on  mooring  line  scope.  Proceeding  from  left  to  right  past  the  knee  of  the  curve 
there  is  only  a  very  slight  decrease  in  cable  scope  required  for  a  large  increase 
in  submergence.  This  is  because  the  mooring  system  is  extremely  taut  (almost 
vertical) .  In  this  regime  the  dominant  force  is  buoyancy ^B,  due  to  submergence 
On  the  left  side  of  the  knee  a  large  increase  in  scope  results  in  a  very  small 
decrease  in  submergence  (i.e. ,  the  buoy  is  approaching  the  freely  floating  case) 
This  regime  is  called  a  slack  line  configuration  where  the  dominant  buoy  force 
is  drag. 

To  determine  the  resultant  buoy  pitch  angle,  the  following  procedures 
would  be  used. 

1 .  For  a  specified  mooring  line  scope  and  current  speed 
determine  the  buoy  submergence,  H,  from  Figure  1-4. 


2.  For  the  same  current  speed  and  for  the  value  of 
submergence  determined  in  Step  1 ,  the  buoy  pitch 
angle  can  be  determined  from  Figure  1-2. 

For  example,  for  a  current  profile  with  a  1 .0  knot  surface  current  and 
for  a  scope  of  2100  feet,  the  submergence,  H,  equals  1.86  feet.  Referring  to 
Figure  1-2  for  a  1.86  foot  submergence  for  a  1.0  knot  current  speed  the  resultant 
buoy  pitch  angle,  e  ,  is  1.93  degrees. 


I 


BUOY  SUBMERGENCE.  H.  IN  FEET 


Figure  1-4.  Buoy  submergence  versus  Mooring  Line  Scope  and 
Current  Profile  Speed 
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DYNAMICS  OF  A  FREE  FLOATING  BUOY 

2 . 1  NATURAL  FREQUENCY  OF  ROLL 

The  natural  frequency  of  roll  (or  pitch) ,  f  ,  for  a  freely  floating  buoy 
is  given  by 

f  -  -L_  J  wjrc 

O  2  TT  f  I 

v 

where 

W  is  the  buoy  weight 

gm  is  the  buoy  metacentric  height 

I  is  the  buoy  mass  moment  of  inertia  about  the  axis 
v 

of  pitch  through  the  center  of  gravity 
The  metacentric  height ,  gm ,  is  given  by 

gm  =  +  gb 

where 

I  is  the  moment  of  inertia  of  the  waterplane  area 

wp 

V  is  the  displaced  volume  of  water 

gb  is  the  distance  between  the  centers  of  buoyancy 

and  gravity 


I 

wp 

I 

wp 

gb 

V 


TT 


4,77 

_2 _ 

4 


4 

) 


25.41  ft4 


9.46  ft 

,R2Ld=  (^f7)2  (50) 
893.5  ft3 
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V 


Therefore , 


—  _  25.41  ft 


gm  = 


+  9.46  ft 


893.5  ft" 
gm  =  9.49  ft 


The  virtual  moment  of  inertia ,  I  ,  is  given  by 


I  =  I  +  I' 
v 

where  I  is  the  buoy  mass  moment  of  inertia  and  l'  is  the  mass  moment  of  inertia 
of  the  entrained  water. 


Using  the  parallel  axis  theorum  and  theorum  of  composite  shapes,  the 
total  buoy  mass  moment  of  inertia,  I,  is  given  by 


where 


i  =  E 

i 


(ii + 


M.  d.2) 

l  l 


— -  th 

I  is  the  principle  moment  for  the  i  component 

til 

Mj  is  the  mass  of  the  i  component 

d.  is  the  distance  between  the  center  of  gravity  of  the 
*  th 

i  component  and  the  center  of  gravity  of  the  buoy 


For  the  aluminum  shell, 


7753  ,  4.77  2  .60^  . 

LA  32.2  1  8  12  ] 

IA  =  72,918  ft  -  lb  -  sec2 
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For  the  sand 


I  = 


M  (-7-  + 
s  4 


12 


I  = 


40000 


32.2 


(  4.77Z  18. 282 

'  1  1  O  / 


16 


12 


I  =  36,359  ft  -  lb  -  sec 
s 


For  the  water  inside  the  buoy, 


D2  L 

,  R  w 

I  -  M  (  ■  +  7—  ) 

w  w  4  12 


=  9416  ,4.772  (26.73  -  18. 28)2  . 

w  32.2  '  16  12  ’ 


I  =  2,156  ft  -  lb  -  sec' 
w 


The  total  moment  of  inertia  is  then 


I  =  I.  +  M.  d  2  +  I  +  M  d  2  +  I  +  M  d  : 
A  AA  s  ss  w  ww 


I  =  226,797  ft  -  lb  -  sec 


The  moment  of  inertia  of  the  entrained  water  is  given  by 

2 


I'  =  M 


D2  L  "  2 

<t  +  -fr  >  +  9b 


,  =  57,169 
1  32.2 


4.77 


50 


16 


12 


)  +  9.41' 


i'  =  531,577  ft  -  lb  -  sec' 


I  =I  +  l' 
v 


I  =  758,374  ft  -  lb  -  sec' 


The  natural  frequency,  f  ,  is  therefore 


(57169  lbs)  (9.49  ft) 
758,374  ft-Lb-sec2 


f  =  0.1346  CPS 
o 


T  =  —  =  7.43  seconds 
o  f 

o 
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2.2 


FREE  FLOATING  BUOY  RIGHTING  MOMENT 


The  righting  moment  acting  on  a  free  floating  buoy  can  be  expressed  in 
terms  of  the  righting  moment  due  to  buoy  tilt,  due  to  the  wave  slope,  and  due  to 
a  combination  of  the  two. 

2.2.1  Righting  Moment  Due  to  Buoy  Tilt 

Righting  Moment  due  to  buoy  tilt  is  expressed  as 

M  =  Wgm  sin  9 

and  is  illustrated  in  Figure  2-1.  The  metacentric  height  gm  is  given  by 

gm  =  bg  +  lim  bm' 

9 — o 

gm  =  bg  +  bm 
I 

—  _  T~  <  wp 

gm  -  bg  +  v 

where  *  is  the  submerged  volume . 

For  small  angles  of  tilt,  the  expression  for  the  restoring  moment 
can  be  linearized  such  that 

M  =  Wgm  9 

2.2.2  Moment  Due  to  Wave  Slope 

The  moment  due  to  wave  slope  as  shown  in  Figure  2-2  is  exactly 
the  same  as  for  buoy  tilt  but  acts  in  the  opposite  rotational  direction. 

That  is , 

M  =  -  Wgm  sin  9 
M  =  -  Wgm  9 

2.2.3  Moment  Due  to  Buoy  Tilt  and  Wave  Slope 

Figure  2-3  illustrates  the  geometry  for  a  combination  of  buoy  tilt 
and  wave  slope.  From  the  previous  expressions  it  follows  that  the  moment 
is  given  by 
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M  =  Wgm  sin  (  9  -  <?) 

M  =  Wgm  (  9  -  9 ) 

2.3  EQUATION  OF  MOTION  FOR  FREE  FLOATING  BUOY 

The  equation  of  motion  for  roll  oscillations  in  waves  is  given  by 

I  V  +  b.  8  +  Wgm  (9-o)  =  o 
v  d 

where  I  is  the  virtual  mass  moment  of  inertial  and  b ,  is  the  linearized 
v  d 

damping  coefficient.  The  wave  slope  9  can  be  expressed  as 

9  s  tan  9  =  a  sin  xt 
where  from  the  theory  of  ocean  waves 


a  =  ILK  =  ^  A 

L  g 

The  equation  of  motion  can  therefore  be  expressed  as 


I  6  +  b  ,6  +  Wgm  6  =  a  Wgm  sin  out 
v  d 

Simplifying  and  substituting  for  constants  we  get 

•  2  2 

(9+2n(9+p^=ap  sin  xt 


2n  = 


p  =  \  =  natural  frequency 

v 


The  steady  state  solution  to  this  equation  will  be 


$  =  c  sin  xt+c  cos  x t 
1  ^ 


Differentiating  we  get 


6  =  Cj  X  cosx-t  -  C2  x  sinx't 

•  •  2  2 
6  =  -Cj  x  sin  xt  -  c2  x  cosxt 
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Substituting  back  into  the  original  equation  and  equating  sine  and  cosine 
terms  we  get: 


,  2  2.  0  2  _  2 
(p  -  a)  )  c^  -  2na>  c  =  a  p 

2ncj  c^  +  (p^  -  a>2 )  ~  0 


2ruo 

C2  2  2  C1 
p  -  (O 

2 


,2  2  .  (2n cj) 

(p  -  oj  )  c1  +  - 


2  2.C1  =  ap 


(P  -  <*  ) 


2.2  2, 

=  a  P  IP  -  ) _ 

1  ,2  2.2  ,  „  .2 

(p  -  oj  )  +  (2na>) 


C2  = 


2ncj 


r  2.2  2. 

a  p  (p  -  oj  ) 


(p"  -  <0  L(p“  -  oj2)2  +  (2na,r  J 


The  solution  can  be  expressed  as 


6=0  sin  (wt  -  a  ) 


Integrating  over  the  length  of  the  submerged  cylinder  we  have 
M  =  1/8  P  CdD  (l*+  i42  )|  9^9 

s  ince 

&  =  0  sin  a>t 

#o 

0  =  00)  cos  ojt 
o 

The  absolute  value  of  the  average  angular  speed,  Q  ,  over  one  cycle  will  be 
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n/2 

JiL  sin  cat 
c o 

o 


Therefore,  since  the  damping  moment  is 

M  =  1/8PCdD  (f4  +  ?4  )  10-1 5 
then  the  linearized  damping  coefficient  is 

b  =  1/8PCdD  (f4  +f4)  |0|0 
For  the  following  parameters 

3 

=  1.99  slugs/ft 


D 


1.5 


D  =4.77  ft 
t  =  15.53  ft 

t2  =  34.47  ft 

(9  =0.2  rad 

o 

The  damping  coefficient  is  computed  to  be 
b  =  333,107  ft-lb-sec 


2.5  STATISTICAL  ROLL  RESPONSE  TO  OCEAN  WAVES  EXCITATION 

The  Response  Amplitude  Operator  (RAO)  is  the  linear  response  to  a  simple 
harmonic  wave  of  unit  amplitude. 


RAO  =  Y  (co  )  =  -|- 


Y  U  )  = 


2  2 
(j _ p 


2  2  2 
g  \/(p  -  co  )  +  (2  ncj  ) 


For  a  wave  height  spectrum,  S  (o>),  the  root  mean  square  of  the  response  can  be 
obtained  from 


R  (  o )  =  Y  (  o  )  S  (  ) 
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The  Pierson-Moskowitz  double  height  spectrum  is  given  by 


m  C 


s2H(w)  =  m_e  “ 

co5 


rl 


*  .  I 


where  V,  is  the  sustained  wind  speed  in  knots, 
k 


Figure  2-5  show  the  buoy  roll  RAO,  Y  (cu)  as  a  function  of  frequency. 
Figure  2-6  shows  the  Pierson  -  Moskowitz  double  height  wave  spectrum  for  a 
wind  speed  of  30  kts  .  The  roll  response  spectral  density  for  a  30  kt.  wind  speed 
is  shown  in  Figure  2-7.  The  area  under  the  curve  in  Figure  2-7  will  give  the  root 
mean  square  response  of  the  buoy  in  roll  for  the  30  knot  wind  speed  spectrum. 

The  process  shown  in  Figures  2-5  through  2-7  was  completed  for  a  range  of  winds 
speeds  and  the  results  are  presented  in  Figure  2-8.  The  three  waves  shown  in 
Figure  2 -8  corres  pond  to 


the  average  of  the  one-third  highest  roll  responses 


the  average  of  all  roll  responses 


0  the  most  probable  (frequent)  response. 
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•  ‘ 


Circular  Frequency,  co,  Radians/Second 

Figure  2-6.  Pierson  -  Moskowitz  Double  Height  Spectrum  for  a  30-knot  Wind  Speed 


Section  3 


U 


DYNAMICS  OF  MOORED  BUOY 


3.1 


RESTORING  MOMENT 


The  restoring  moment  as  illustrated  in  Figure  3-1  is  given  by 

+ 

~2 


L.  +  H  \ 

M  =  (B  +  AB)  gm  (ft-<p)  +  T  L  9-  Dn  -2-* -  -  L  -  T  L 

3  y  w  °  \  z  w  /  x  w 


where , 

gm  =  ( 


Ln  -  L  J  B  +  (l 
B  w/  \  o 


+  #-  -  L  )aB 


w/ 


I 


B  +  AB 


+  wp 


(l„  -  L  ')b  +  /l  +—  -L  )  AB  I  _ 
\  B  w  /  \  o  2  w/  wp  7 


gm 


B  +  \3 


B+  AB 


M  = 


H 

+  o  -  L 


(LB-Lw)3  +  (Lo  +  ^-Lw)AB+TI 


wp 


{9-9) 


+  T  L  9  -  D, 


y  w 


/  L.  +  H  \ 

(  — -  -  L  )  -T  L 

\  2  w/  x  w 


M  = 


wp 


(VLW)  B  +  (Lo+  f-Lw)AB+Tl 

(LB-gB  +  (Lo+r  -R>+^ 


wp 


+  T  L  d  -  Dn  feJliL  -  L  )  -  T  L 
yw  o  \  2  w  /  xw 


The  angle  ^  is  the  sum  of  the  steady  state  and  time  varying  angles, 
ft0  and  tf(t) 


9  =  ft  +  (9 (t)  =  ft  +  ©sin (cut -cr) 
o  o 


•  • 

9=  ft  (t)  =©cucos  (cut  -cr) 

V  =  ft  (t)  =  -Qu?  sin  (cut  -cr) 
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ul 

4 


4 

t 


jB  +  V 


M  =  L  -  L  |  B  +  IL  +  —  -  L  AB+7I  +  T  L  \0 
B  w j  \o  2  w)  1  wp  y  wl  c 


o  2  wl 


wp  y  w  o 


+  L  -L  B  +  L  +  2  -L  AB  +7 1  +  T  L  &(t) 


wp  y  w 


-  VLw 


1  B  +/l  +  f  -l  Iab  +7l  o 

I  \  o  2  w/  wp  ■ 


>  (-** 

°\  2 


-  L  -  T  L 

W  /  X  w 


The  steady  state  pitch  angle,  SQ,  is  given  by 

0  [  D0  (.b.+  B..\ 

“o  =  sin0Q  =  V  2  j 

BLb  +  AB  (  0  *Z~  ] 


BLb  +  AB 


Lr’L 
B  w  j 


)  B  +  (lo  +  f  -L  )AB  +7l  +T  L  I - 2- 

1  1  °  w/  wp  y  w|blb  +ab(lq+ 


o'  2 


+  [(3  +  A  B)  gm  +T  L  1(9  (t)  -  (B  +  A  B)  gm  p 
1  y  w 


/  L  +  H 

-Do(-V- 


r  /  lq  +  h  \ 

=  (l  -  L  )b+(l  +“p*) AB  +  7l  - - 1 - 

L  jLb(lb-lJ +^B(Lo+r 


+  (B  +  A3)  gm  +  T  L  fi(t)  -  (B+\B)  gm  ? 
I  y  w1 

/L  +  H\ 


r  - 

is 
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EQUATION  OF  MOTION 


The  equation  of  motion  is  given  by 


I  9  +  9  +  [(B+AB)  gm  +T  L  1  6  (t)  =  (B 

v  d  L  y  wJ 


9  =  a  sin  cut 


a  = 


JlH  =  a 

L  g 


3  +  AB  =  W  +  Ty 

6  + 


o ,  • 

A)  i  t) 


(W+  T  )  gm+T  L 
y  y  w 


i 


(W  +  T 


tfft)  = 


V 


8+aii+bd  =c  sin  cut 
b. 


a  = 


I 

v 


b  = 


c  = 


(W  +  T  )  gm  +  T  L 
y  y  w 

.  —  “ 

v 

(W  +  T  )  gm  2 

- ^ A 

I  g 


The  steady  state  solution  is  given  by 
H-  sin  cut  +  cos<ut 

Differentiating  we  get 

0  =  coscut  -  c^cu  sincut 

2  2 
d  =  c,  u  sincut  -  c„  -u  coscut 


+  AB)  gm 


)  gm 

_  a  sincut 


Substituting  back  into  the  original  equation. 


-CjO>  sino)t  -  cos  u>  t  +  ac^o  cos  cat  -ac^osinc^t 

+  be  ^  sin  cot  +bc2C0st->t  =  csinojt 
Equating  coefficients  of  sine  and  cosine  terms , 

-c^oj  -  ac^u  +  bc^  =c 
2  . 

- c 2<^>  +  ac^  +  bc2  =  o 


(b  —  oj  )  c^  -  au>  —  c 


a<uj  c1  +  (b-  oj  )  =  o 


-aoj 

C2  “  2  C1 

b  -oj 


(b-  oj  )Cj  + 


b-  u 


2  C1  =  C 


2  2  2 
(b  -  a)  )  +  (a  u) ) 

b  -  oj2 


c^  =  c 


c,  = 


C  (b  -  w  ) 


(b  -  a;2)  +  (a  u  j  ' 


sin(wt  -  a) 
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The  Pierson-Moskowitz  double  wave  height  spectrum  is  given  by 

,4 


-9.7  x  10 
135  V?  .4 


SM  =  e  *  cj 

ur 


The  response  spectral  density  is  given  by 


R(w)  =  Y  U  S  U 


R(oj)  = 


(W  +  T  )2  gm2  135 

y 


^(W  +  Tv)gm+TvLw  A  2  /b,' 


2  _  2| 

ajg  Iv 


-9.7  x  10 
—5 — g— 

e  V  w 


•A2  /b,\ 

.  _ y  w  z\  I _d  I 

v  i  “  \i  ) 

v  v 


-/ 


R  =  /  R(cj)  dcj 

*0 


In  order  to  evaluate  these  functions  it  is  first  necessary  to  relate 
several  mechanical  parameters  to  the  buoy  submergence  term,  H,  since  the 
buoy  is  no  longer  considered  freely  floating.  The  parameters  which  must  be 
computed  are  the  new  virtual  mass  moment  of  inertia,  the  new  linearized 
damping  coefficient,  and  the  effective  metacentric  height. 

Calculate  new  virtual  mass  moment  of  inertia 

I  =  I  +1' 

v  new  new  new 

I  =  I  ,  .  =  226,797  ft-lb-sec2 
new  old 

I*  I'  , ,  +  Ul'  +  AM  d., 

new  =  oid  \  w  o>  / 
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=  t»  +  a 

new  old  g 


+  (io  +  2  j 


I’  =  531,577  ft-lb-sec2  +  35.56  H  1 . 42  +  +  Iz4 . 47  + 

new  1 2  \ 


Calculate  new  linearized  damping  coefficient,  b^ 

b  =  7  -|r-  P  CL  D  (  L  4  +  (L  -  L  +  H)4) 
d  4  7T  T)  \  w  o  w  ’ 


15.53  ft 


50  ft 


4 . 77  ft 


1.99 


lb- sec 


0.2  rad 


bd  =  ^  ^  (1.99)  (1.5)  (4.77)  (15.53)4  +  (50  -  15.53  +  H)4 


bd  =  0.227  (58,168  +  (34.47  +  H) 


Calculate  new  Metacentric  Height 

m  -  (1b  1  LA  B  +  (Ib*  f  -  Lwu  B  j  !wP 

B  +  A  B  TTJZ  (I^  +h) 


=  (24.99-  15. 53)  (57169)  +  (50  +  ^  -15.53hB  + 


25.41 


57169  +  PB 


540819  +  (34.47  +t)^B 


57169  +  A  B 


7H4.77 

4 


17.87(50  +H) 


If  the  following  mooring  configuration  is  assumed 

•  1.0  knot  surface  current  profile 

•  2080  ft  scope 

•  2  ft  buoy  submergence,  H 

•  1.9  degree  steady  state  pitch,  Qq 

then  the  following  parameters  can  be  calculated  to  be 


AB  =  Ty  =  2290  lbs 

I  =  848,047  ft-lb-sec2 

v 

b,  =  414,072  ft-lb-sec 

d 

gm  =  1 0 . 4  9  ft 

Using  the  above  values,  the  Roll  Response  Amplitude  Operator, 

Y(uj)  can  be  determined  and  is  shown  in  Figure  3-2.  The  Pierson  Moskowitz 
double  wave  height  spectrum  for  a  30  knot  wind  speed  is  shown  in  Figure  3-3 
Figure  3-4  shows  the  roll  response  spectral  density  to  the  30  knot  wave 
spectrum.  Figure  3-5  shows  the  statistical  mean  roll  response  as  a  function 
of  wind  speed. 


009 


Circular  Frequency,  ^ ,  in  Radians/Second 

Figure  3-3.  Pierson  -  Moskowitz  Double  Wave  Height  Spectrum  for  a  Typical  30  Knot  Wind  Speed  as  a 
Function  of  Frequency 


Wave  Height  Spectrum 


Moored  Buoy  Mean  Roll  Response  as  a  Function  of  Wind  Speed 


